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Abstract 

In this paper, we show numerically that a symmetric Earth matter density profile can, in principle, be recon- 
structed from a single baseline energy spectrum up to a certain precision. For the numerical evaluations in the 
high dimensional parameter space we use a genetic algorithm. 
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1. Introduction 

In order to obtain knowledge about the Earth's 
interior, current measurements are chiefly based 
on seismic wave propagation through the Earth 
(see, e.g., Refs. The methods used in geo- 

physics allow a quite precise reconstruction of the 
seismic wave velocity profile. However, for the 
determination of the matter density profile many 
assumptions have to be made about the equation 
of state, relating the velocity and matter den- 
sity profiles. This process involves several uncer- 
tainties, since additional information about the 
Earth's interior is quite rare (see, e.g., Refs. |^,^). 
In neutrino physics, a method similar to X-ray to- 
mography has been suggested, namely neutrino 
absorption tomography |5Hl3||. It requires many 
different baselines, since neutrino absorption is 
not sensitive to the arrangement of the matter 
structure along one baseline. Moreover, the cross 
section rises with neutrino energy, which means 
that high energetic neutrinos are required. 

Neutrino oscillations in the Earth have been 
extensively studied in general and it has been 
suggested that neutrino oscillations in matter 
can severely alter the energy spectrum of a neu- 
trino beam through the Earth for long baselines 
|l4|-^. For quite realistic calculations of the 
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transition probabilities the Preliminary Reference 
Earth Model (PREM) or similar matter density 
profiles have been used |]3T| , ^ . In addition, sev- 
eral approximations of the matter density profile, 
such as a step function (see, e.g., Refs. |]3^ , ^) o r 
the first terms of a Fourier series expansion [|35| , 
turned out to supply rather good results. How- 
ever, the inverse problem, i.e., what we can learn 
about the Earth's structure from neutrino oscil- 
lations, is also considered to be an interesting 
problem. In Refs. (3^,^^ the approach of an in- 
verse scattering problem was suggested. However, 
transition probabilities of a single baseline energy 
spectrum turned out not to be sufficient for deter- 
mining the matter density profile uniquely, i.e., 
additional information on the relative phase of 
the output state vector was required. In this pa- 
per, we ignore this information and investigate 
what we can learn about the Earth's matter den- 
sity profile from a single neutrino baseline en- 
ergy spectrum using two flavor neutrino oscilla- 
tions in matter. We show that we can, in prin- 
ciple, reconstruct the (symmetric) matter den- 
sity profile, up to a certain precision without any 
additional assumptions about the profile itself. 
The (large-scale) symmetry of the matter density 
profile seems to be a reasonable assumption by 
the results from geophysics. Furthermore, since 
two flavor neutrino oscillations cannot distinguish 
time reverted matter density profiles (see, e.g., 
Refs. |38| , |39| ), omitting the symmetry condition 
would a priori mean that reconstructed, asym- 
metric profiles are not unique. 
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2. The model 

For the numerical evaluation we use a (non- 
deterministic) genetic algorithm |^,^. It as- 
sumes an initial generation of random matter den- 
sity profiles and calculates their respective energy 
spectra. The energy spectra are compared to a 
(realistic or measured) reference spectrum by an 
error (fitness) function. In this case, we choose 
minimization of a x^-function 

B 

X' = V2[(x,)-x,]+2a:,log-^, (1) 

where B is the number of energy bins, Xi the (in- 
teger) number of events in the ith bin, and (xi) 
the mean number of events in the ith bin. Non- 
deterministic evolution of the matter density pro- 
files over time then creates an optimal or subop- 
timal matter density profile with an energy spec- 
trum similar to the reference spectrum. Since the 
whole process is completely nondeterministic, the 
algorithm may find any symmetric matter density 
profile which fits the energy spectrum with the 
required precision. Thus, it also indicates if we 
are really able to obtain a unique solution for the 
given energy spectrum. The strength of applying 
a nondeterministic method is to reduce the ex- 
ponential calculation effort to a polynomial one, 
where the quality of the obtained result can be 
determined by the statistical x^-analysis. This 
means that we will find samples of the (high di- 
mensional) parameter space within a (high di- 
mensional) ncr-contour. 

Here we use the neutrino appearance chan- 
nel of a single neutrino^ baseline with an en- 
ergy range from 2.5 GeV to 20 GeV (an energy 
range which may be produced by a neutrino fac- 
tory), in order to cover the Mikheyev-Smirnov- 
Wolfenstein (MSW) resonance regions for reason- 
able data p^"|45|] . Since we are in principle inter- 
ested in the reconstruction of the Earth's mat- 
ter density profile, we ignore information about 
the beam energy distribution and the cross sec- 

^ Since for Am^ > antineutrinos do not show resonant 
behavior in matter, they are much less suitable for a re- 
construction of the matter density profile. Thus, we will 
restrict this discussion to neutrinos. 



tions, i.e., we use the transition probabilities di- 
rectly. However, it is useful to have a statis- 
tical estimate of the significance of the results. 
Thus, we assume a total number of events I to 
be folded with the transition probabilities equally 
spread over the whole energy spectrum at an 
equal cross section. Furthermore, we use a two 
flavor neutrino oscillation scenario with a sym- 
metric PREM profile ||3|,||. For the propaga- 
tion of a neutrino state vector along the base- 
line L through the Earth, we divide the (sym- 
metric) matter density profile into 2N layers of 
equidistant length AL = ALk = L/{2N) of con- 
stant matter density pk, where fc = 1, 2, . . . , 2N, 
and we use an evolution operator in flavor ba- 
sis in each individual layer Uf{ALk, Pk) (see, 
e.g., Refs. [ ^6p^ ). The total evolution opera- 
tor U f (L) is then the time-ordered product of 
the ones in the individual layers, i.e., Uf{L) = 
Uf{AL2N, P2n) ■ . . Uf{AL2, p2)Uf{ALi,pi). Fi- 
nally, the transition probabilities are given by the 
absolute values squared of the elements of the to- 
tal evolution operator, e.g., the transition prob- 
ability for Vf3 is Pap{L) = \ {vi3\Uf{L)\Va)\- 

Note that, because of the symmetry, N is the 
number of independent parameters in this prob- 
lem. The electron density ne^k is related to the 
matter density pk by 

rie.fc = ^, fc = l,2,...,2iV, (2) 

where Yk is the average number of electrons per 
nucleon in the fcth layer. In the Earth, Yk — 
^Earth = 1/2, k = l,2,...,2iV, and ttin is 
the nucleon mass. For the baseline we choose 
L = 11736 km, corresponding to a nadir an- 
gle of about 0.4, in order to be able to see the 
Earth's mantle as well as core in the energy spec- 
tra. Furthermore, for the oscillation parameters 
we use Am2 = Aml^ = 3.2 • lO'^eV^ @ and 
sin 26' = sin 26*13 = 0.1 [|9|,|o|, which are applica- 
ble for a mass hierarchy Arn^i ^ Am^2- 

3. Results and analysis 

We now present the results of some genetic al- 
gorithm trial runs. For the number of steps N 
in the matter density profile we choose = 4, 
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iV = 9, and = 14. We also compute the ref- 
erence energy spectra with the PREM profile di- 
vided into 2N layers of constant matter density 
for = to be equal to zero at the minimum. 
In almost every example, we use a total number of 
events / = 100000 to be folded with the transition 
probabilities, corresponding to a very optimistic 
guess of a very large neutrino factory, as well as 
a number of energy bins i? = 31. However, we 
will also investigate the precision for a varying 
number of events and energy bins. 

The best fit results of the genetic algorithm are 
shown in Fig. |l| for several numbers of steps TV of 
the matter density profile. Note that the refer- 
ence energy spectrum has two peaks correspond- 
ing to the mantle at high energy and the core as 
well as interference effects at low energy, respec- 
tively. All results in this figure are within a la 
(68%) confidence region, which indicates that we 
cannot measure the matter density profile pre- 
ciser without additional assumptions. Therefore, 
the precision of the measurement is limited for 
this type of experiment. Nevertheless, every trial 
run of the algorithm converges to a matter density 
profile clearly showing the mantle-core edge with 
high quality. It could not be seen in a Fourier ex- 
pansion approach, since many Fourier coefficients 
had to be determined to obtain a high resolution 
of this edge ||. 

The fact that the initial profiles were com- 
pletely created at random, indicates that a sin- 
gle neutrino baseline energy spectrum really de- 
termines the arrangement and matter densities 
of the large-scale structure of the Earth's matter 
density profile. Especially, in comparison to X- 
ray or neutrino tomography, it can distinguish the 
order of the different layers of matter because of 
the non-commuting operators used for the propa- 
gation through the matter density profile. The al- 
gorithm may have found any matter density pro- 
file matching the reference profile, such as one 
with mantle and core exchanged. Nevertheless, it 
turned out that it converged in all cases against 
the well-known mantle-core-mantle structure of 
the Earth. However, the non-commuting opera- 
tors make the analytical inversion of the energy 
spectrum quite complicated (c/., Refs. p6| , ^ ). 

As already mentioned above, the resolution of 



the matter density profile seems to be bounded. 
In order to investigate this, let us take a closer 
look at the parameter space. Since it is impossi- 
ble to visualize the contours of equal confidence 
levels of a high dimensional problem, we plotted 
in Fig. H for = 14 some representatives of the 
parameter space within the Icr, 2(t, and 3(t con- 
tours. Again, the mantle-core edge can be easily 
resolved at the Id-level, but in many cases not for 
lower precision. It turns out that this effect grows 
with an increasing number of steps in the matter 
density profile. In addition. Figs. and ^ indi- 
cate that the matter density profile can be quite 
well reconstructed for a small number of steps N, 
but for a higher number of steps N small fiuctu- 
ations enter the reconstructed profile. As can be 
seen in Fig. ^ we especially observe for iV = 14 
small fluctuations in the mantle and core regions, 
indicating bounds on the spatial resolution. In 
App. we show analytically with a perturba- 
tion theoretical method what is intuitively clear: 
structures of small amplitude on a length scale A 
much smaller than the oscillation length in mat- 
ter L°f^ttcr = 47r£;/ATO2 = ATTE/[£,{n^)/\m'^] can- 
not be resolved by neutrino oscillation in matter. 
Here 
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cos 26* 



sin-" 29 



(3) 



is determined by the mean matter density n,,. 
Since ^ > |sin20|, the upper bound A ~ 
-^matter < 47r£'/(sin 26'Am^) oc E indicates that 
the higher the resonance energy of the resonance 
peak is, the lower the spatial resolution becomes. 
For this reason, the resolution of the core, cor- 
responding to the low energy peak, is probably 
higher (c/.. Figs. |l| and ||). However, for a real- 
istic neutrino factory one has suppressions in low 
energetic events, which means that this effect is 
probably compensated by statistics js^. 

In order to investigate the dependence on the 
total number of events / and energy bins B, we 
show in Fig. ^ some matter density profile ex- 
amples close to the Icr-contour for different num- 
bers of / and B. One can see that the precision 
increases with the number of events. In almost 
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Figure 1. The best fits of the genetic algorithm trial runs for iV = 4. = 9, and = 14 steps in the 
Earth's matter density profile. On the left-hand side the best fit matter density profiles (solid lines and 
boxes) are plotted with the reference profiles (dashed lines and stars). The values of N and are given 
above each plot. On the right-hand side the energy bin spectra of the best fits (light bars) are shown 
together with the spectra of the reference (dark bars). 
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Figure 2. Representatives of the = 14-dimcnsional parameter space for la (first row: < 15.9), 
2cr (second row: ^ 24.0), and 3a (third row: < 33.2) confidence levels. The respective x^'s of 
the samples are given above each plot. The sample matter density profiles are plotted with solid lines, 
whereas the reference profile is plotted with dashed lines. 



all cases, we can clearly resolve the mantle-core 
edge, but for too few events I <^ 5000 the result 
does not have a physical meaning anymore. One 
can also see that the precision increases with the 
number of energy bins B. In addition, for B = 5 
the matter density profile again looses its physi- 
cal meaning. Note that, since we did not integrate 
over the energy of each bin, but used the mean 
energy, there is a statistical error inc;reasing with 
decreasing number of energy bins. However, one 
may expect a natural boundary of the evaluation 



for B < N, because A'' represents the number of 
unknowns and B the number of equations. This 
equivalence between the number of equations and 
the number of energy bins would be destroyed for 
an integration over the energy within each energy 
bin. 

4. Summary and conclusions 

In summary, we have demonstrated that we 
can, in principle, reconstruct the symmetric 
Earth matter density profile using a single neu- 
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Figure 3. Matter density profiles for different numbers of total events / (left) and different numbers of 
energy bins B (right). The figures show matter density profiles close to the Icr-contour (x^ ~ 10.42) for 
= 9 steps. The following matter density curves arc plotted: In the left plot {B — 31): Reference (thick 
soUd); / = 1000 (thin soHd); / = 5000 (dotted); / = 10000 (dashed); / = 100000 (dashcd-dottcd). In the 
right plot (/ = 100000): Reference (thick solid); B = b (thin solid); B = % (dotted); B = 16 (dashed); 
B = 31 (dashed-dotted). 



trino baseline energy spectrum without addi- 
tional assumptions. However, it turned out that 
the precision is limited in this method, especially 
since we made quite optimistic assumptions about 
the total number of detectable events / and the 
energy spectrum of the source, as well as we ig- 
nored any energy smearing. 

Wc conclude that neutrino physics cannot be 
used as the only source for obtaining information 
about the Earth's interior, since the precision is 
presently lower than in geophysics, and the res- 
olution has a natural lower bound, the oscilla- 
tion length in matter. However, geophysics can 
only access the velocity profile of seismic waves 
directly. Additional assumptions about the equa- 
tion of state have to be made in order to access 
the matter density profile. Thus, neutrino oscil- 
lations in matter could help to test and verify the 
equation of state. Similarly, the parameter space 
in neutrino reconstruction tomography could be 
shrinked by using knowledge from geophysics. 
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A. Spatial resolution of neutrino oscilla- 
tions in matter 

In this appendix, we will show in the case of two 
neutrino flavors that we cannot resolve structures 
much smaller than the oscillation length in mat- 
ter. For that we state the differential equation 
describing the time evolution of flavor states in 
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matter: 



= Ho Ufit) + 7^i„t(ne) Vf{t) (4) 



^ sin 261 ^ cos 26* 



and 
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Let us introduce a perturbation 77 (t) on a shorter 
timescale than a slowly varying matter density 
profile fieit), i.e., 



ne{t) = ?"ie(0 + ^7(0 



(7) 



with 77 (i) making many oscillations while rie(t) is 
approximately constant. This will also result in 
a modification of the amplitudes 



(8) 



where i'f{t) solves the unperturbed equations for 
'^int(^e)- Thus, wc Can split Eq. by applying 
Eqs. and (^) and we obtain 



(9) 
(10) 



Here we have neglected second order corrections 
of the order i]{t)xit). We know that the os- 
cillation of the solution of the first equation is 
determined by the oscillation length in matter 
L-^tter = ^7:E/Am^ = AnE/[^{n,)Am% where 
the f was defined in Eq. and refers here to 
the constant mean matter density rie.^. For the 
second equation, we may assume a quick periodic 
oscillation ri{t) = Acos{ujt) with the amplitude A 
and the frequency w on a timescale much shorter 



*For the approximation of a matter density profile by its 
mean, see Ref. |47[. 



than the one of the first equation. Hence, we 
can in the second equation take £'/(t) — I'f and 
he(t) ~ n approximately to be constant. 
We can thus rewrite Eq. as 



noXfit) + nintine)xfit) 



(11) 



This is a differential system of equations simi- 
lar to the one for two coupled oscillators, where 
the one representing the electron flavor is driven 
by an external force corresponding to the oscilla- 
tory fluctuation in the matter density proflle. It 
is solved by quite lengthy, purely oscillatory ex- 
pressions for Xe{t) and Xnit)- 

One can show that the amplitudes of 
the oscillations reduce to zero when uj ^ 
Amy{2E) = 27r/L°^l,,,,. Since ^ = 2^/A, 
this means that we cannot resolve structures in 
the matter density profile with a width A <^ 



T ° 



Note that the perturbation method used 



breaks down for too large amplitudes or too slow 
oscillations. 
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